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Abstract 

Field distributions generated by static QQ and QQQ sources are calculated analyti- 
cally in the framework of Gaussian (bilocal) approximation of Field Correlator Method. 
Special attention is paid to the QQQ system and asymmetric configurations are also 
studied. At large quark separations fields form distinct stringlike shape. In both cases 
QQ and QQQ the string consists mainly of longitudinal color electric field. Transverse 
color electric component contribution is shown to be less then 3%. Baryon string has 
an y-like shape with a deep well at the string junction position. Field distributions 
for quark-diquark and for three quarks on one line are considered. The interaction 
potential for quarks forming an equilateral triangle is calculated. The material of the 
paper is illustrated by 24 3D colored pictures. 

1 Introduction 

Field distributions inside the string connecting static QQ sources have been measured re- 
peatedly on the lattice using both connected [1-3] and disconnected ||, |5| probes. Similar 
measurements were done later also for Abelian projected configurations ||. Analytic calcula- 
tions for the disconnected probe made in in the framework of the Gaussian approximation 
to the Field Correlator Method (FCM) [EL [J have revealed a clear string-like structure of 
the same type as was found on the lattice. The connected probe allows to measure the sign 
and modulus of each field component separately. 

A comparison of lattice data [p]] on magnitude and direction of fields in the QQ system 
with analytic predictions of FCM was done in [|2| , demonstrating a remarkable agreement in 
all distributions. In particular, the measured decrease with distance of longitudinal electric 
field from the string axis ("the string profile") was remarkably well described by contribution 
of the lowest, (bilocal) correlator [[J. It should be noted that the input of FCM is the the 
field correlator, defined by scalar formfactors D and D x ||. The lattice measurements [10] 
yield for the latter the exponents with the slope T g « 0.2fm. The dominance of the bilocal 
correlator (sometimes called the Gaussian Stochastic Model (GSM) of the QCD vacuum) 
was verified recently on the lattice by the precision measurement of static potential for QQ 
Wilson loop in different SU(3) representations [|11|]. Analysis of data JTl| made in [0] has 
demonstrated that GSM contribution around 99% to the static QQ potential is consistent 
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with the data. These results give an additional stimulus to the analytic calculation of field 
distributions using lowest bilocal correlator. 

In this paper we calculate QQ and QQQ field distributions in bilocal approximation of 
Field Correlator Method using connecting probes. We have found that the string consists 
mainly of longitudinal color electric field. For three quarks forming an equilateral triangle the 
string has an F-like shape with a deep well at string junction position. We consider also field 
distributions for quark-diquark and for three quarks at one line cases. Besides, we calculate 
three-quark-interaction potential for quarks forming an equilateral triangle. The well around 
the string-junction position provides a decrease of slope of the interaction potential at small 
quark separations. 

Some results of this paper were presented earlier in [13]. 

The paper contains 5 sections. In section 2 the QQ longitudinal and transverse field 
components distributions are derived analytically and plotted in 3D pictures for different 
quark separations. Total field distributions including perturbative one gluon exchange are 
also presented. In section 3 we derive analytically three quark field distribution and plot 
it for different quark configurations. For three quarks forming an equilateral triangle we 
plot also total field including perturbative part. In section 4 baryon interaction potential is 
calculated. In section 5 the results are summarized and their possible physical consequences 
are discussed. 



2 Field distributions in meson 

We shall study the field distributions using gauge invariant construction p^(x) consisting 
of probing plaquette P /Jil/ (x) and Wilson loop W connected with parallel transpoters $ (see 
Fig.l). In what follows this construction is referred to as a connected probe. By definition, 

f v (W^x )^(x ,x)(P, u (x))](^)j(x,x )} 
PA X ) = — " 1 Jy^ !. I 1 ) 

where 



W = ^ t rPexp(^£^ a ^), (2) 

(P»„(x)r p = (Pexptga 2 F; u (x)t a )^ (3) 
^(x,y) = (Pex pi g f ' ApTdz^, (4) 

Wff (x ) in (1) denotes an ordered exponential along contour C of Wilson loop without the 
point Xq = (0,0,0,0). A rectangular contour C of size RxT lies in the plane (1 4) (since 
quarks are static). The surface S is bounded by the contour and has coordinates x' = 
(x'i, x f 2 , x f 3 , x' 4 ), where < x[ < R, x' 2 = x' 3 = 0, — T/2 < x' A < T/2. The probing plaquette 
Pfiu{x) of size a x a is oriented in the plane (// v). It is placed at a point x = (xi,X2, x^, X4), 
where X\ is the probe coordinate along the string axis, x 2 is the distance from the probe to 
the string axis , x 3 = rr 4 = 0. In the small-plaquette-size limit expanding (3) in powers of 
a 2 one obtains that the connected probe is propotional to F(x): 

pAx) = iga — ^ + 0{a ). (5) 
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Thus using the connected probe allows to measure QQ color field components not disturbed 
by probing plaquette in the small a limit. We shall calculate the connected probe (1) 
in bilocal approximation of FCM expanding Wilson loop in F pv and keeping only bilocal 
correlators contribution. Let us write Wilson loop in terms of using the nonabelian 
Stokes theorem: ^ 

W = —exp(ig jf F pu (x, z )da fiu (x)), (6) 

where 

Fftu(,X, Zq) = $(Zq, x)F flv (x)$(x, Z ). (7) 

Averaging Wilson loop over vacuum fields in bilocal approximation one has 

(W) «exp(-^^d^(x)d^ ff (x0^tr(F^(x)$(x,xOF pCT (xO$(x^2:))^ (8) 

In bilocal approximation we arrived at a double surface integral. The connected probe 
surface S p consists of Wilson loop surface S and plaquette surface Sp, S p = S + Sp, which 
gives for the double integral in (8) the sum of three terms, 



. , +2 / / + / / . (9) 

S J S J S ^ Sp J S p ^ S p 

Using (9) and (8) for the connected probe (1) we realize that the first term in (9) cancels 
with the denominator in (1) and third term of (9) disappears in the small a limit, yielding 

, exp(-|J gp J Sp ) _ 
exp(-£/s/ s ) 



exp( 



, \\ I ')-! = - / / +0(a 4 ), (10) 

S j Sp Z J Sp J Sp J S J Sp 



where 



ISJSp * 



= -a 2 j^da pa (x')—tr(F pcT (x')§(x',x)F pv (x)§(x,x')). (11) 
Let us define now an averaged (colorless) field strength at point x in bilocal approximation 

as 

(F pu (x)) QQ ee -^l oc {x) = 



j da pa {x')^tr{F pa {x')®{x\x)F^{x)<S>{x,x')). (12) 



c 2 

da pa (x')- T 

is 

In FCM the following parametrization of bilocal correlators by scalar formfactors D and Di 
is suggested (see second ref. in |§): 

^tr(F pa (x')$(x',x)F pi/ (xMx,x')) = (5 Pfl 5 al/ - 5 pu 5 ap )(D(h 2 ) + D^h 2 )^ 

3D (h 2 ) 

+ (h p h p 5 ua - h p hj up - h p h u S pa + h v hj pp ) — — — = D paiflu (h), (13) 
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where h = x — x' . 

In lattice measurements JTH] both formfactors were found to be exponential beyond x = 
0.2fm with a slope T g ~ 0.2fm: 

D(h 2 ) = D(0) exp(-(i\h\), D x {h 2 ) = D x (0) exp(-/i|/i|), (14) 



D 1 (0)~±D(0), //-lGeV; T 9 = -. 

We will use the form of D, £)i (14) in the whole region of x as it was done in M. 

The values (p cr) = (l 4) on the rhs of (12) are determined by the Wilson loop orientation, 
and therefore the combinations of D and D\ entering (12) are: 

DuAh) = S u (D(h 2 ) + D 1 (h 2 )) + {hxhi + hl5 u ) ^ 2 \ (15) 



Du,ik(h) = {hkh^n - hih 4 5 k i) — ^ — , (16) 



dD^h 2 

where i,k — 1,2,3. One can see that (15) gives color electric field components (E(x))qq, 
and (16) the color magnetic ones (B(x))qq. 

Let us show that the only nonzero field components are (Ei(x))qq and (E2(x))qq. Really, 
in (15) DuM h ) = for h=x 3 - x' 3 = 0, so that in (12) (E 3 (x)) qQ = (F u (x)) qQ = 0. 
Du } ik{h) in (16) is antisymmetric in h A = —x' 4 , so after integration in (12) over J s da pa (x') = 
J dx\ J_ T / 2 dx' 4 one obtains (B(x))qq = 0. 

Let us calculate (Ei(x))qq and (E 2 (x))qq. 



(Mx))qq= L~dx[ I^Jx'^Dih^+D^ + ihl + hl)^^) = (17) 



r R r T/2 ( 

\ dx\ / dx\ D{h 2 ) + D x {h 2 ) + (h\ + h\ 

JO J-T/2 \ 



At T -> oo 

rxi poo , 

(Ei(x))%q = D{0) / dhi / dh A exp{-^hj + h 2 2 + hj 

Jxi—R J — oo 



i xi— R 

fXl 



2D(0) r dxJx 2 + xlK 1 (fiJx 2 + xl), (18) 

Jx x -R V V 

where K\ is McDonald function. 

D(0) 
3 

£>(0) 



J dh A (xi exp(— p^xf + x 2 + hi) — (x x — R) exp(— n\J Jx\ — R) 2 + x\ + h\ 
x x sjx\ + x\K x (\Lsjx\ + x 2 2 ) - ( Xl - R)\]{X! - Rf + X^K^fly/ (x x - R) 2 + X$f) . 



(19) 



{E 2 {x)) Q q = dx 1 dx 4 h x h 2 — 7^3 — = / dh x dh A - 

W Jo J-T/2 Oh 2 Jxx-R J -00 



t/2 dh 2 Jxx-R J-00 2 dhi 



D(0)x 2 
3 

D(0)x 2 



J dhi ^exp(— n\j x\ + x\ + hi) — exp(— ji\J (x\ — R) 2 + x\ + h 



v y x? + x\K x {^x\ + x|) - y/( Xl - Rf + xlK 1 {^{x l - Rf + x|)J (20) 
Let us proceed with some evaluations. At R T s 

/■-R/2 

E?(R/2,0) =4D(0)J o dxxKi{x/T g ) = 2irD(0)T g 2 (21) 

At quark position 

(0, 0) = 2D(0) / dxxKAx/T g ) = vrD(0)T o 2 = a, (22) 

JO 

where cr=0.9GeV/fm is the QQ string tension. The same value of a one has in the bilocal 
approximation of FCM from law of Wilson loop. 
The fields E 2 and E? 1 are maximal at x ~ T 9 , 



£ 2 (0,T 9 ) = ^(r„0) = ^Zjl^l) = °-^D(0)T* (23) 



We conclude that E™ x /E™ x « 0.38/12 = 3.2%. 

Perturbative interaction in the leading a s order is defined by the one gluon exchange 
between quarks which leads to the Coulomb-type field contribution, 

E Coul = C F a s - ^ , (24) 

where ri is a vector from the quark to the probe, and r 2 is from the antiquark to the probe. 
For fundamental quarks and QQ in SU(3) singlet the Casimir operator Cp = trt a t a = |. 
The parameter e = |a s = 0.295 is defined using the lattice mesurements of QQ potential 
fitted by Cornell potential Vcom = — 5 + crR (see review [14] and references therein). 

In Figs. 2-5 (a)Q we plot (Ei(xi,x 2 ))qq distributions for quark separations R = T g , 5T g , 10T 9 
30T g . One can observe in the pictures how the string of a characteristic shape between quark 
and antiquark is forming starting from R = 5T g . On Figs. 2-5 (b) the total field distributions 
with perturbative one gluon exchange included are plotted for the same separations, 

(E&) 2 = (^ Coul + (E^qq) 2 + (£ 2 Coul ) 2 . (25) 

One can see that the one gluon exchange dominates at separations smaller then T g = 0.2fm. 

In Fig. 6 (E 2 (xi, x 2 ))qq distributions for R = T g and 30T 9 are presented. One can prove 
that E 2 does not create by itself the string. The magnitude of E 2 at distances of the order 
of T g around the quark and antiquark positions is less then 3% of £j nax and is decreasing 
fast with distance from Q (Q). 

In Fig. 7 are presented the transverse crossections of distributions at fixed x\ = R/2 
plotted in Figs.2-5(a). For R = !0T g and 30T 9 the string profiles practically coincide and 
at R > 10T g the string acquires its saturation form. The field magnitude in the middle of 
string is E sat = E™* = 1.8GeV/fm= 2a. The string width is 5x SSjt = 2.2T g (by definition of 
5x sat , Ef(R/2, <5x sat /2) = l/2(£ sat ) 2 ). 



L You can get quality ps figures from [http://hcron.itcp.ru/~kuziricnko/figurcs.tar.gz 
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In Fig. 8 a longitudinal string crossection at fixed x 2 = for R = 30T g is shown. The field 
is increasing rapidly in the range from — 3T g to 3T g and then is forming the long plateau 
Ei = E sat . In the same figure the total field with perturbative one gluon exchange included 
is shown. It is decreasing monotonically at x\ > 0, flattens at x\ ~ 5T g and is increasing at 
25T 9 < xi < 30T g ; in the whole range of < Xi < R sign of derivative does not change. 

In Fig. 9 the field distributions (Ei(xi,X2))qq are presented for R = T g and 30T g . Just 

as E 2 , Ei 1 does also not create by itself the string. Magnitude of E® 1 at distances ~ T g 
around quark and antiquark positions is less then 3% E sa,t and is decreasing fast off the 
quark (antiquark). One can distinguish two symmetrical distributions around quark and 
antiquark. At R = T g they are superposed and the distribution has a maximum at zero. At 
R = 30T g a region in which the field is absent appears between them. These distributions are 
antisymmetrical under transformations Xi — > —xi and (xi — R) — > — (xi — R) respectively. 

3 Field distributions in baryons 

In this section we study field distributions generated by three static quarks in different 
configurations. Let us first consider three quarks forming an equilateral triangle. The contour 
of baryon Wilson loop W^ 3Q ^ consists of three contours C r , T = A,B, C, formed of quark 
trajectories, and the string junction trajectory (see Fig. 10). The string junction position is 
determined by the minimal area condition. In Fig. 10 three contour lopes are intersecting 
at angles forming a Mercedes star. The gauge invariance of the loop is provided by the 
antisymmetric tensors: 

= le^ 7 $«,(C A )^(C B )^,(C c )e^' (26) 

The probing plaquette P^ v (x) is attached by parallel transporters to the contour Ca at a 
point X . X -y — X ^ = 0, ^3 = R. For the entire connected probe construction the coordinate 
x 2 = 0. 

In the expression of averaged field in the bilocal FCM approximation (12) we are now to 
sum over three surfaces A, B, C : 




Here Sp is the surface corresponding to the contour Cr, T = A,B,C; the sides of the 
rectangular surface have an extention R and T. Introducing in the plane (xi,xs) unit vectors 
n A = (0, 1), n B = (^y, —\), n c = (^f^, — §), we integrate in (27) over the surfaces: 

da T pc {x')F p(T {x') = dal 4 (x')F i4 (x') = < E t (l'n r , x' 4 )dl'dx' 4 (28) 

with the result 

(F^(x)) 3Q = 5>[ f R dl' [ T/2 dx' 4 D l4 ^(l'n r - x,^ - x 4 ). (29) 

r J J-T/2 

In what follows we shall only consider the D contribution which, as was demonstrated in 
previous section , gives the string form. From (29), (15), (14) one obtains 

(E 3 (xi,x 3 )) 3Q = D(0) / dl' / dx' A {e-^ + ^- 1 '^- 

JO J-T/2 
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/P / x\ n / n N /" T/2 A I <^> -wj (xi-^l') 2 +(x3+l'/2) 2 +x'i 

{Ei(xi,X3))aQ = D(0) / dZ / cte 4 (— e V 2 

JO J-T/2 2 



-T/2 

_ ^_ e ~^\l {xi + ^V?+(x 3 +V/2) 2 +x'i^ ^ 

As one can see from (30), (31), the baryon string is a superposition of the meson string 
(Ei(x 3 , Xi))qq given in (18) and two such meson strings obtained by rotation over angles ^ 



and f : 



(E 3 ( Xl ,x 3 )) 3Q = (E 1 (x 3 ,x 1 ))° Q - ^i(4^'i))go " \(M4,^)qq, (32) 
(Ei{x u x 3 ))3Q = ^(E^x'J)^ - ^(E.ixlx'i))^ (33) 



where 



= -2 Xl ~ "y^ 3 ' x 3 = ~Y x i - 2 X s; (34) 

1 a/3 „ a/3 1 
= -2 Xl + ~y X3 ' x 3 = — 2~ Xl ~ 2 X3 ' 

In Figs. 11-13 (a) and Fig. 14 the squared field distribution is plotted 

(E(xi, x 3 ))Iq = (E l {x 1 ,x 3 )) 2 3Q + (E 2 ( Xl ,x 3 )) 2 3Q (35) 

for i? = T g ,5T g , lOT g and 30T r A characteristic feature of all distributions is a deep well 
around the string junction position. At the junction position itself the field vanishes, since 
for the symmetry reason at this point a preferred direction is absent. Outside the well the 
baryon string becomes a sum of three meson strings going from the quarks to the string 
junction. Let us note that field distribution in the region around string junction will remain 
intact even for different distances from quarks to string junction were all the "outer" meson 
strings saturated. 

The Coulomb field of one gluon exchange for three quarks is 

= + % + <*» 

where is distance from the i-th quark to the probe. Factor ^ appears due to contraction 
of antisymmetric tensors in Wilson loop (26): 

1 1 Cp 

—^afi^aar^fi a ^ 1T = ~Q^ T ^ O ~ ^~ ' (3^) 

where Cf = | is Casimir operator; Cfol s = 0.295 (see explanation after (24)). The total 
field 

Eg = (E>30 + Eg ul (38) 

is given in Figs.ll-13(b) and Fig. 15. In Fig. 11(b) the total field is at least two orders of 
magnitude greater than the nonperturbative part (Fig. 11 (a)) in entire region considered. In 
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Fig. 12(b) we still see a half of string, going out of a Coulomb spike. In Fig. 13(b) and Fig. (15) 
the Coulomb field only changes the shape of the adjacent piece of the string and does not 
disturb the rest of it. 

In Fig. 16 are shown the longitudinal crossections of strings along x 3 axis at x\ — 
depicted in Figs.ll-13(a), 14. At R > 5T g the shape of the well does not depend on 
quark separation. At R > 10T 9 the string becomes saturated; its longitudinal crossection 
at < x 3 < 6T g rises from to (_E sat ) 2 and then levels up. In the negative x 3 region 
< \x 3 \ < 1.25T 9 the crossection grows from to 0.25GeV 2 /fm 2 and decreases rapidly at 
1.25T 9 < |x 3 | < 5T g . Let us define a well radius R wcU as E 2 (0, R wcU ) = (E sat ) 2 /2. Then from 
the Fig. 16 we find R wcll = 1.75T g . 

Let us consider now a quark-diquark configuration when two quarks are close to each 
other and far from the third one. We take them to form an isosceles triangle with a base 
much shorter than its lateral sides. The string is formed along the minimal path, consistsing 
of three segments connecting quarks with string junction and intersecting at angles like 
a Mercedes star. Let us denote Rqq the length of two short segments and Rq that of the 
long one. In Eq.(29) to calculate the field distributions one can now integrate over dl 1 in the 
range from to Rq for V = A and from to -Rqq for T = B,C. Adding to Eqs.(32),(33) 
the indices corresponding to integration ranges, one has 

{E 3 {x u x 3 )) Q _ QQ = (E 1 (x 3 ,x 1 ))^ - \{E^, x[))^ - 

-\{E,{xlx'[)) D Q ^, (39) 

(E 1 (x 1 ,x 3 )) Q . QQ = ^{E^x'Xp* - ^(E 1 (4,x'Xq QQ - (40) 
In Fig. 17 we plot distributions 

(E(xi, x 3 )) 2 Q _ QQ = {E 1 (x u x 3 ))l_ QQ + (E 3 ( Xl ,x 3 )) 2 Q _ QQ (41) 

for Rq = SOTg and two values of -Rqq, (a): Rqq = 0.5T g , (b): Rqq = 3T g . In case (a) one 
observes a mesonic string (cf. Fig. 5 (a)). In case (b) one can see a baryonic string with the 
well in the string junction area. How this transition happens? In accordance with (39), (40) 
the field strength at x\ = x 3 = is determined by the difference of the field values of mesonic 
strings of length Rq Rqq at this point. The field strength of the mesonic string (18) E (R) 
at the origin for the string of length R is 

E Q (R) = — / K^xdx. (42) 

7T JO 

This function increases linearly from zero at small R and saturates to the asymptotic value 
-E'asymp = o" at -R > 4T g . The dependence Eq(R) is shown in Fig. 18. 

Let us define a radius of transition of baryonic string into mesonic one -Rbar^mes as 
E^(R har ^ mes ) = El symp /2. From Fig.18 we find i? b ar^mcs = 1-5T S . At R QQ > R har ^ mcs the 
quark-diquark string has a characteristic baryonic form with the well at the string junction 
(the case of Fig. 17(b)) and at Rqq < -Rbar->mcs it turns into a mesonic one (the case of 
Fig-17(a)). 

Let us now consider the case when quarks are placed along one line at distances R\ and _R 2 
between them. The string junction position, as in previous cases, is defined by the minimal 
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string length condition. For given quark positions the string length is minimal when the 
string junction position coincides with that of a middle quark (placed between two others) 
and is equal to the sum of distances from this quark to others. 

The contour of Wilson loop is in the plane (14). Proceeding just as in the beginning of 
section, we introduce in the plane (1 3) two unit vectors ni = (— 1,0), n 2 = (1,0) and get 

(E 3 ( Xl ,x 3 )) 3Ql = 0, (43) 

Ri rT/2 

dl' I dx' 4 D 1414 (-l' - x 1} —x 3 , x' 4 ) + 

J-T/2 

rK 2 i rT/2 / f 

+ dl' dx' 4 D 141 M' - xi, -x 3 , x 4 ) = 

JO J-T/2 

= -( El {- Xl ,x 3 )) D Q f + (E^x^f. (44) 

Distributions (E 1 (xi,x 3 )) 3 q 1 for Ri — R 2 — 15T 9 and R\ = 10T g ,R 2 = 20T 9 are shown in 
Figs. 19 (a),(b). Since both mesonic strings have saturated profiles, the field at the string 
junction position x 1 = is exactly zero. 



(E l (x l ,x 3 )) 3Ql = - J 



4 Quark interaction potential in baryon 

In this section we consider static quarks placed in vertices of an equilateral triangle at 
distance R from the string junction. 

As the Wilson loop is a Green function of three static quarks, it gives the quark interaction 
potential in the baryon, 

V^\R) = - hm — \n{W ( - 3Q \R, T)). (45) 

T— >oo X 1 

Using bilocal approximation we get an expression of potential as a sum of surface integrals 
of bilocal correlators over Wilson loop surface. We integrate over the surfaces just as in 
calculation of averaged baryon field. 

-(W^\R,T))= 1 - £ i fda; u (x)da b pa (x')x 

Z a,b=A,B,C Jlja ,7Lf > 

x ^-tr{F^(x)$(x, x')F pa (x')^(x', x)F^(x)) = 

3 rR rR rT/2 r T/2 

= - V nfn b k / dldl' / / dx 4 dx' 4 D l4k4 (ln A - l'n b ,x 4 - x' 4 ). (46) 
2 b=A B C Jo Jo J-T/2 J-T/2 

In the last equation of (46) we have used symmetry of the Wilson loop contour - 

E =3 E ■ (47) 

a,b=A,B,C a=A; b=A,B,C 

As in the previous section we shall only consider the contribution to the potential from 
the formf actor D. Then 

3D(0) f* f* JIJf , 

-T/2 J-T/2 



Let us call the first exponential in (48) a "diagonal" as it is obtained at a = b = A in (46) 
and second one — "nondiagonal" as it is obtained at a = A, b = B,C in (46). Integrating 
in (48) one obtains 

y(3Q) (jR) = f R dUK^l) - -1(2 - (piR) 2 K 2 (piR))- 

7T JO jJ, 6 




+ + V/2y Kl {^-^ + (l + l>/2)*)}, (49) 

where we have used the normalization of D(0) (22). 

First two terms of (49) result from the diagonal and the third one from the nondiagonal 
exponentials of (48). The first term yields linear potential with a slope 3er at distances 
R ^> T g , the second is a small correction to the first one (at large R it equals VlaTg/ir). The 
third term is increasing and significantly contributing to the potential growth till R ~ 3T g 
and then flattens. This term is due to the well in the middle of baryonic string. 

Potential behaviour V( 3Q \R) is shown in Fig.20(a). Tangent slope cr'(R) = %(R) in- 
creases from zero to 3<r in the range of < R < QT g . Beyond this range one can consider 
baryonic string as a sum of three mesonic ones with the radius of well is subtracted (cf. Fig. 8 
and Fig. 16). 

It Fig.20(b) the total potential V^\r) is presented with the one gluon exchange in- 
cluded, 



VS?\R) = £ - + V™(R) = -^f^ + V^{R). (50) 

rij = Ry^3 is quark separation, summation is performed over three pairs, CVa s = 0.295 (see 
explanations after (24), (36)). In the figure one can see a change of the potential slope as a 
consequence of the well in the middle of baryonic string. 

Let us compare results of our calculations with lattice measurements of baryonic potential 
(review [14] and refs. therein). In the latter the potential for quarks forming an equilateral 
triangle was fitted by a Coulomb plus linear (Cornell) potential in the range of 0.055fm< 
R < 0.71fm. As a result the slope of the linear potential was found to be 2.6a which is equal 
to the a'(3.5T g ). Thus the effect of the potential decrease due to the well in the middle of 
baryon string is consistent with the effect obtained on the lattice. 



= £ - + V^\R) = ^ + V^(R). 



5 Summary 

We have calculated connected-probe field distributions for QQ and QQQ cases using the 
lowest (Gaussian) field correlator. Since the Gaussian correlator contributes 99% to the 
static QQ potential [11,12], one expects that higher correlators would not change our picture 
significantly. 

The connected-probe analysis allows to distinguish directions of components of color 
electric and color magnetic field due to probe orientation. One can see that the string 
mainly consists of the longitudinal color electric field. The transverse electric component 
contributes less than 3% of the longitudinal one. The magnitude of longitudinal component is 
given by the formfactor D and the transverse one by D\. The D 1 contribution to longitudinal 
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component is less than 3%. At separations R greater then lOT g the QQ string profile becomes 
saturated and does not change with increase of R. The saturated string width is equal to 
2.2T g . Our QQ results are in agreement with earlier calculations in [1-3]. The bulk of the 
string fields is also in accordance with disconnected-probe analyses 0, ||. 

The baryon string picture is obtained in our paper for the first time. A deep well in the 
electric field distribution appears around the string-junction position, where electric field 
vanishes because of symmetry arguments. The well has a radius of 1.75T g demonstrating 
strongly suppressed fields in the middle of the heavy baryon. For quark-diquark configuration 
we have defined the radius of baryonic to mesonic string transition, which was found to be 
equal to 1.5T g . For three quarks on one line the field at the middle quark position was found 
to vanish when quark separations are large. Since the Wilson loop for the QQQ configuration 
used to generate interaction is also applicable for light quarks fl5| , physical consequences of 
this well can be in principle observable both for light and heavy hadrons. One consequence 
can be illustrated by Fig. 20, where the nonperturbative part of the potential grows very 
slowly at small R, so that the asymptotic slope is obtained only at very large distances. 
Therefore an effective slope for ground-state hadrons can be some 10-20% smaller, the fact 
which is in agreement with relativistic quark model of baryons [[HJ and with recent lattice 
calculations of static QQQ potential fH| . 

One should note that the vanishing of fields at the string junction holds for the connected- 
probe analysis of each separate field component, and may not be true for field fluctuations 
at this point, measured by the disconnected probe. 

Financial support of the grants 00-02-17836 and 00-15-96786 is gratefully acknowledged. 
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Figure 1: A connected probe for QQ. 

Figure 2: (E 1 (x 1 ,x 2 ))qq (a) and (Eq^(xi, x 2 )) 2 (b) distributions in GeV 2 /fm 2 . Quark 
separation R = T g . Both x\ and x 2 are measured in T g . Q and Q positions are marked with 
points. 

Figure 3: (E 1 (x 1 ,x 2 )) 2 q q (a) and (E t Q t Q (x 1 , x 2 )f (b) distributions in GeV 2 /fm 2 . R = 5T g . 
Both xi and £2 are measured in T g . Q and Q positions are marked with points and vertical 
lines. 

Figure 4: The same as in Fig. 3 but for R = 10T g . 
Figure 5: The same as in Fig. 3 but for R = 30T g . 

Figure 6: (E 2 (x u x 2 )) q q distribution in GeV/fm for R = T g (a) and R = 30T g (b). Both 
xi and x 2 are measured in T g . Q and Q positions are marked with points and vertical lines. 

Figure 7: QQ string transverse crossection profile {Ei(R/2, x 2 ))qq in GeV 2 /fm 2 . x 2 is 
measured in T g . Dotted, dashed, long dashed and solid lines refer to R = T g , 5T g , lOT g and 
30T g . At R larger then lOT g the string becomes saturated; R = lOT g and R = 30T g profiles 
practically coincide. 

Figure 8: (E^xi, 0)) 2 QQ (solid line) and (Eg£(xi, 0)) 2 (dashed line) in GeV 2 /fm 2 . Xl is 
measured in T g . R = 30T g . 

Figure 9: D\ contribution to (E\{x\, x 2 ))qq distribution in GeV/fm for R = T g (a) and 
R = 30T g (b). Both x 1 and x 2 are measured in T g . Q and Q positions are marked with 
points and vertical lines. 

Figure 10: A connected probe for QQQ. 

Figure 11: (E(xi, x 3 ))\q (a) and (E%q(xi, x 3 )) 2 (b) distributions in GeV 2 /fm 2 for quarks 
forming an equilateral triangle. Quark separation from the string junction R = T g . Both x\ 
and X3 are measured in T g . Quark positions are marked with points. 

Figure 12: (E(x u x 3 )) 2 Q (a) and (Ef^(x 1: x 3 )) 2 (b) distributions in GeV 2 /fm 2 . R=5T g . 
Both X \ and £3 are measured in T g . Quark positions are marked with points and vertical 
lines. 

Figure 13: The same as in Fig. 12 but for R = 10T g . 
Figure 14: The same as in Fig. 12(a) but for R = 30T g . 
Figure 15: The same as in Fig. 12(b) but for R = 30T g . 

Figure 16: QQQ string profile (E(0 : x 2 ))Iq in GeV 2 /fm 2 . x 3 is measured in T g . Dotted, 
dashed, long dashed and solid lines refer to R — T g , 5T g , 10T g and 30T g respectively. 

Figure 17: (E(x\, x 3 ))q_qq distributions in GeV 2 /fm 2 for quark-diquark configuration 
with Rq = 30T 9 ,Rqq = 0.5T g (a) and Rq = 30T 9 ,Rqq = 3T g (b). Both x\ and x 3 are 
measured in T g . Quark positions are marked with points and vertical lines. 

Figure 18: E 2 (R) = ((£i(0, 0))g o ) 2 (i?) distribution in GeV 2 /fm 2 as a function of quark 
separation R measured in T g . 

Figure 19: (E\(x\,xz))\q 1 distribution in GeV 2 /fm 2 for R\ = R 2 = 15T g (a) and R± = 
10T g ,R 2 = 20T g (b). Both X \ and x 3 are measured in T g . Quark positions are marked with 
points and vertical lines. 

Figure 20: Three-quark-interaction potential V^®\R) (a) and the total potential with 
perturbative one-gluon-exchange included V^\R) (b) as functions of distance R from 
quarks to string junction measured in T g . 
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